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First-principles study of ferroelectricity and pressure-induced 
phase transitions in HgTiOa 

Alexander I. Lebede\0 

Physics Department, Moscow State University, 119991 Moscow, Russia 
(Dated: February 25, 2013) 

Ground-state structure is found and pressure-induced phase transitions up to 210 kbar are studied 
in mercury titanate from first principles within the density functional theory. It is established that 
the _R3c structure experimentally observed in HgTiOa is metastable at ambient pressure. With 
increasing the hydrostatic pressure, the ground-state structure changes following the _R3 — >■ i?3c — >■ 
Pbnm sequence. It is shown that the appearance of ferroelectricity in HgTiOs at P = is associated 
with an unstable phonon mode. Optical and elastic properties of different phases of mercury titanate 
are calculated. The quasiparticle band gap calculated in the GW approximation [Eg = 2.43 eV) 
agrees with experimental data better than the value obtained in the LDA approximation (1.49 eV). 
Analysis of the thermodynamic stability explains why the synthesis of mercury titanate is possible 
only at high pressures. 

PACS numbers: 77.84.Dy, 61.50.Ks, 64.70.Kb, 77.80.-e 
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I. INTRODUCTION 

The strain engineering has already become an impor- 
tant technological approach which enables to enhance the 
properties of many electronic materials. In ferroelectrics, 
due to the large strain-polarization coupling, the strain 
effects are especially important. For example, in ferro- 
electric thin films and superlattices, the biaxial strain 
can be used to finely tune the ferroelectric, dielectric, 
and piezoelectric properties of these materials fi"— 

The influence of hydrostatic pressure on displacive 
phase transitions was first analyzed by Samara et al.,'^ 
who explained why an increase in pressure usually de- 
creases the temperature of phase transitions associated 
with soft optical phonons at the center of the Brillouin 
zone and increases the temperature of phase transitions 
associated with soft phonons at the boundary of the Bril- 
louin zone. Recently, a very different behavior of the fer- 
roelectric properties was discovered in PbTiOa at high 
pressures; the observed enhancement of the ferroelectric 
instability with increasing pressure was explained by the 
original electronic mechanism of ferroelectricity.™ The en- 
hancement of ferroelectricity at high pressures was pre- 
dicted for many oxides with the perovskite structure i^ 
These interesting findings necessitate further studies of 
the pressure effects on the ferroelectric properties. 

Despite limited experimental data on mercury ti- 
tanate, this material exhibits interesting but contradic- 
tory ferroelectric properties. Mercury titanate HgTiOs 
can be prepared from HgO and Ti02 at pressures of BO- 
GS kbar,M The obtained crystals have a rhombohedrally 
distorted perovskite structure. The observation of second 
harmonic generation (SHG) in HgTiOs at 300 ¥J enabled 
to propose that mercury titanate is non-centrosymmetric 
and its space group is RZc. However, because of limited 
accuracy, the atomic coordinates in Ref. 'TI were deter- 
mined only for centrosymmetric structure i?3c. Sub- 
sequent studies of dielectric properties of mercury ti- 
tanate^iS did not found sharp dielectric anomalies: a 



strongly asymmetric broad peak with a maximum dielec- 
tric constant of '^800 at about 220 K and a noticeable 
hysteresis in the heating-cooling cycle as well as weak 
narrow peak at about 515 K were observed. At 300 K 
no dielectric hysteresis loops were observed in electric 
fields up to 10^ V/m..?'^ Scanning calorimetry revealed 
weak anomalies in the 420-480 K temperature rangef^iS 
but their temperatures differed from the temperatures of 
maximums in dielectric measurements. 

X-ray diffraction studies of HgTiOa under hydrostatic 
pressure^ revealed non-monotonic behavior of the do24 
interplanar distance and of the (104)-(110) doublet split- 
ting at P « 20 kbar, which was explained by a phase 
transition from the rhombohedral to the cubic phase. 
Studies of the electronic structure^° of rhombohedral and 
cubic modifications of HgTiOa using the full-potential 
linearized augmented plane wave (FP-LAPW) method 
showed that the rhombohedral Ric phase is a direct- 
gap semiconductor with the band gap energy of ~ 1.6 eV, 
whereas the cubic phase is a metal. 

To resolve the inconsistency of the ferroelectric prop- 
erties of HgTiOa and predict other properties of this 
material, first-principles calculations of the ground-state 
structure, phonon spectra, optical and elastic proper- 
ties of mercury titanate at hydrostatic pressures up to 
210 kbar were performed in this paper. 



II. CALCULATION DETAILS 

The calculations were performed within the first- 
principles density-functional theory (DFT) using pseu- 
dopotentials and a plane-wave expansion of wave func- 
tions, as implemented in ABINIT software. ^^ The lo- 
cal density approximation (LDA) for the exchange- 
correlation functional was used. Optimized separable 
nonlocal pseudopotentials^ were constructed using the 
OPIUM program,!'^ The parameters used for the construc- 
tion of pseudopotentials are given in Table U (for meaning 



TABLE I. Parameters used for construction of pseudopotentials. All parameters are in Hartree atomic units except for the Vioc 
energy which is in Ry. 



Atom 



Configuration 



Td 



Qp 



Id 



Viae 



Hg 


sd^oes^p" 


1.78 


2.00 


1.78 


7.37 


7.07 


7.37 


Ti 


3s'3p«3d''4s° 


1.48 


1.72 


1.84 


7.07 


7.07 


7.07 





2s2 2/3d» 


1.40 


1.55 


1.40 


7.07 


7.57 


7.07 



0.01 



1.41 



2.65 



of the parameters see Ref. |12I ): for Ti atom the local po- 
tential correction^- was added. 

The pseudopotentials for Ti and O atoms were non- 
relativistic (they have been already tested and used in 
Ref. \b), the pseudopotential for Hg was constructed 
using scalar-relativistic generation scheme. The plane- 
wave energy cut-off used in the calculations was 30 Ha 
(816 eV). Integration over the Brillouin zone was per- 
formed on 8x8x8 Monkhorst-Pack mesh for the dielec- 
tric phases and 12x12x12 mesh for the metallic phases 
with one formula unit in the unit cell. For larger unit 
cells, equivalent k-point density was used. The total en- 
ergy was converged to less than 10~^° Ha. The structure 
relaxation was stopped when the Hellmann-Feynman 
forces were befow 5 x 10~^ Ha/Bohr (0.25 meV/A). 

The convergence tests with a number of irreducible 
k-points increased by 2-4 times have shown that the 
changes in the calculated total energies did not exceed 
0.05 meV for the dielectric phases and 0.35 meV for the 
metallic phases. The changes in the phonon frequencies 
were less than 0.01 cm^^ for hard modes and less than 
1 cm""'^ for soft modes. An increase in the energy cut-off 
also had a little effect on the calculated data: the change 
in the relative energies of different phases was less than 
0.15 meV when the cut-off energy was increased to 35 Ha. 

To test the quality of pseudopotential for Hg, the calcu- 
lations for orthorhombic and rhombohedral polymorphs 
of HgO were performed. Among them the orthorhombic 
modification (the montroydite mineral) had the lowest 
total energy. The calculated lattice parameters of these 
phases (a = 3.4663 A, 6 = 6.6253 A, c = 5.3013 A for 
the orthorhombic phase and a = 3.5092 A, c = 8.5417 A 
for the rhombohedral phase) were in reasonable agree- 
ment with experiment's (a = 3.5215 A, 6 = 6.6074 A, 
c = 5.5254 A; a = 3.577 A, c = 8.681 A). Phonon spec- 
tra were calculated with the same interpolation scheme 
as was used for other titanates of group H elements in 
the periodic tableJ^ 

The quasiparticle band gap in HgTiOs was calculated 
using the so-called one-shot GW approximationfiiiS 
The Kohn-Sham wave functions and energies calculated 
within DFT-LDA were used as a zeroth-order approxi- 
mation. The dielectric matrix eGG'(q, i^) was computed 
for 4x4x4 q-mesh from the independent-particle polariz- 
ability matrix PQQ,(q, w) calculated for 4285 reciprocal- 
lattice vectors G(G'), 42 occupied and 158 unoccupied 
bands. The contribution of higher-lying bands was taken 
into account using the approach proposed in Ref. ,19l The 




FIG. 1. Phonon dispersion curves for HgTiOs in the cubic 
PmSm phase. The labels near the curves denote the symme- 
try of unstable modes. The absence of LO-TO-splitting at 
the r point is caused by metallic character of this phase. 



dynamic screening was described using the Godby-Needs 
plasmon-pole model. The components of wave functions 
with kinetic energy below 24 Ha were used in these calcu- 
lations. The energy correction to the DFT-LDA solution 
was computed as diagonal matrix elements of S — E^c 
operator, where E = GW is the self-energy operator, 
E^c is the exchange-correlation energy operator, G is the 
Green's function, and W = e~^v is the screened Coulomb 
interaction. In the calculations of E, the components of 
wave functions with kinetic energy below 24 Ha for both 
exchange and correlation parts of E were used. 



III. RESULTS 

A. Ground-state structure at P = 

The calculated phonon spectrum of HgTiOs in the per- 
ovskite cubic phase (space group PmSm) is shown in 
Fig. [1] It is seen that two types of instability appear 
simultaneously in this phase: the stronger one associ- 
ated with the deformation and rotation of the oxygen oc- 
tahedra (the r25-X3-M3-r25-i?25-.^3 branch) and the 
weaker one associated with the ferroelectric (antiferro- 
electric at the boundary of the Brillouin zone) instability 
(the Ti5-X^-M^-Ti5 branch). The absence of LO-TO- 
splitting at the r point is caused by the metallic band 



TABLE II. Energies and volumes per one formula unit for 
different distorted phases of HgTiOa at P = 0. The energy of 
the cubic phase is taken as the energy reference. The phase 
with a minimum specific energy and minimum specific volume 
are denoted by bold values. 



Unstable 


Space 


Energy, 


Volume, 


mode 


group 


meV 


A^* 


— 


Pm3m 





57.573 


X3 


P4:2/m7nc 


-88 


57.275 


ri5 


R3m 


-94 


58.923 


Tis 


PAmm 


-122 


59.444 


Fas 


P4m2 


-139 


57.088 


ri5, Fas 


Amm2 


-151 


59.749 


X5 


Pmma 


-202 


57.916 


Xs 


Cmcm 


-306 


57.867 


r25 


R32 


-467 


56.956 


R25 


lA/mcm 


-778 


56.188 


Ms 


PA/mbm 


-809 


56.195 


R25 + M3 


Pbnm 


-936 


55.853 


R25 


Imma 


-940 


56.099 


R25 


R3c 


-974 


56.336 


A2u 


R3c 


-982 


56.632 


— 


R3 


1059 


60.140 



structure of cubic HgTiOa. 

To determine the structure of the ground state, the en- 
ergies of different distorted phases originating from the 
cubic structure were calculated for each of the above- 
mentioned unstable modes taking into account their de- 
generacy. As follows from Table |TT1 the R3c phase has the 
lowest energy among these phases. This phase is derived 
from the Pm3m structure by out-of-phase rotations of 
the oxygen octahedra around three cubic axes as a re- 
sult of condensation of triply degenerate R25 mode at 
the boundary of the Brillouin zone {a~ a~ a~ tilt system 
according to Glazer's notation). The energy of this phase 
is even lower than that of the Pbnm phase, in contrast 
to other titanates of group II elements in the periodic 
tableJ^ It should be noted that the overlapping of the 
conduction and valence bands disappears as the struc- 
tural distortions become larger, and all phases with the 
energy lower than —300 meV are semiconductors. 

The ferroelectric instability specific for the reference 
Pm3m structure of HgTiOa still exists in the R3c phase. 
The calculations show that two unstable modes with ^211 
and Eu symmetry and frequencies of 135i and 21i cm~^ 
are observed at the F point in the phonon spectrum of 
this phase. Among the corresponding ferroelectrically 
distorted phases, the i?3c phase has the lowest energy. 
All phonon frequencies at the center and at A, D, and Z 
points on the boundary of the Brillouin zone in this phase 
are positive, the matrix of elastic moduli (Sec. IIII C[) is 
positive-definite, and so the i?3c phase is the ground- 
state structure. The calculated lattice parameters and 



atomic coordinates in R3c and i?3c phases are given in 
Table Hill For comparison, the experimental data for the 
-R3c phased are also included into this table. It is seen 
that the calculated data for the i?3c phase agrees well 
with the experimental data. The calculated interatomic 
distances for the _R3c phase and calculated mean dis- 
tances for the i?3c phase agree with interatomic distances 
determined from X-ray diffraction^ too (Table HV)) . 

Along with the phases originating from the perovskite 
structure, other possible structures, in particular, the 
ilmenite one characteristic for titanates of group II 
elements— MgTiOa, ZnTiOg, and CdTiOg^— should be 
tested. The calculation showed that at ambient pressure 
(P — 0) the ilmenite structure of HgTiOs with R3 space 
group has the lowest energy among all considered phases 
(Table |TTl . The fact that R3c or R3c structures are ob- 
served in X-ray diffraction enables to suppose that these 
phases are metastable. This metastability is evidently 
associated with large difference between i?3c (i?3c) and 
R3 structures, both in the lattice parameter and in the 
rhombohedral angle ( Table Hill . The phase transition be- 
tween so different structures should be of the first order, 
for which a wide region of metastability is characteris- 
tic. The reason why all investigated samples had the 
metastable R3c (i?3c) structure is that the samples have 
been prepared at 60-65 kbar, at which the _R3c phase 
is thermodynamically stable (Sec. IIII D|) . The energies of 
two more possible hexagonal phases, one with the 2H 
BaNiOs structure and other with 6H BaTiOa structure 
(both have PG^/mmc space group), are, respectively, by 
269 meV and 73 meV higher than the energy of the cubic 
Pm3m phase. 



B. Origin of ferroelectric instability 

We now discuss the ferroelectric properties and nature 
of the ferroelectric phase transition in HgTiOa. As the 
changes in the Hg-0 bond lengths accompanying the fer- 
roelectric phase transition do not exceed 0.1 A and the 
difference in energy of the R3c and R3c phases is only 
8.1 meV, the Curie temperature in HgTiOa cannot be 
too high and will not exceed 300 K. Therefore we can 
associate the Curie temperature with the temperature of 
the first maximum in the dielectric constant, which was 
observed at 220 K.*'^ This interpretation is confirmed by 
the absence of dielectric hysteresis loops at 300 K. The 
fact that SHG signal was observed in Ref. at 300 K 
can be explained by the presence of defects; the cause 
of the defect formation in HgTiOs will be discussed in 
SecHlEl 

Calculation of the dynamical matrix for the R3c phase 
shows that the zz-element of the on-site force constant 
for Hg atoms is small but positive (0.0263 Ha/Bohr^), 
and so these atoms cannot be considered as off-center 
ions. The analysis of the eigenvector of the ferroelectric 
^211 mode in this phase shows that the displacements of 
Hg atoms in this mode is 22 times smaller than that of 



TABLE III. The lattice parameters and atomic coordinates in HgTiOa phases with R3c, R3c, and R3 space groups at P = 0, 
and in the Pbnm phase at 141 kbar. 



Phase 


a, A 


a, deg. 


Atom 


Position 


X 


y 


z 


7?3c 


5.4984 


58.4093 


Hg 


2a 


0.24904 


0.24904 


0.24904 








Ti 


2a 


-0.00333 


-0.00333 


-0.00333 








O 


6b 


0.66598 


-0.15240 


0.25846 


7?3c 


5.4881 


58.4252 


Hg 


2a 


0.25000 


0.25000 


0.25000 


(calc.) 






Ti 


2b 


0.00000 


0.00000 


0.00000 











6e 


0.65983 


-0.15983 


0.25000 


i?3c 


5.4959 


58.59 


Hg 


2a 


0.25 


0.25 


0.25 


(exp.)" 






Ti 


2b 


0.0 


0.0 


0.0 











6e 


0.665 


-0.165 


0.25 


Ri 


5.8304 


53.9320 


Hg 


2c 


0.36869 


0.36869 


0.36869 








Ti 


2c 


0.84974 


0.84974 


0.84974 











6f 


0.55966 


-0.03220 


0.19275 


Pbnm 


5.2678 (a) 


— 


Hg 


4c 


-0.00445 


0.03190 


0.25000 




5.2983 (b) 




Ti 


4b 


0.50000 


0.00000 


0.00000 




7.5501 (c) 




O 


4c 


0.08502 


0.47264 


0.25000 











4d 


0.69594 


0.30216 


0.04431 



" The experimental data taken from Ref. J^ were recalculated for rhombohedral setting. 



TABLE IV. Interatomic distances in HgTiOa phases with i?3c 
and i?3c space groups. 



Pair of 


Distance, A 




Number 


atoms 


This work 




Ref. 7 


of bonds 




Ric 


Rlc 




Hg-0 


2.198 


2.195 


2.20(4) 


3 


Hg-O 


2.698, 2.888 


2.786 


2.77(4) 


3+3 


Hg-0 


3.172 


3.162 


— 


3 


Ti-0 


1.906, 2.064 


1.977 


1.96(4) 


3+3 



Ti atoms. This means that collective displacements of 
titanium. atom.s against oxygen atom.s are responsible for 
the ferroelectric instability. Weak ferroelectric activity 
of Hg atoms is confirmed by small values of their Born 
effective charges: Z*^ = Z*^^ = 3.20, Z*^ = 2.42; they 
slightly exceed the nominal charge of the ion. For Ti 
atoms the corresponding values are Z*^ — Z*y — 7.85, 
Zl. = 7.92. 

The calculated static dielectric constant at K in the 
i?3c phase is almost isotropic {sxx = 97, £zz = 101).— Its 
value is compatible with a maximum dielectric constant 
of ~800 observed in the experiment at 220 K4i^ The cal- 
culated spontaneous polarization in the RUc phase turns 
out unexpectedly large, Ps = 0.37 C/m^. Apparently, 
this is a result of large effective charge of the A2u mode 
in the paraelectric Ric phase (Z*jj = 12.66). 



C. Optical and elastic properties of HgTiOa 

The electronic structure calculations performed in this 
work within the DFT-LDA approach confirm earlier re- 
sult^*' obtained in the generalized gradient approximation 
(GGA) that the Ric phase is a direct-gap semiconduc- 
tor with the band gap of 1.6 eV and the band extrema 
located at the T point. In our calculations, the band gap 
of HgTiOs at P = is ^^da ^ ^ 49 gy ^^^ j^-g pressure 

coefficient is dE^^^/dP = +0.44 meV/kbar. However, 
both these results disagree with the experimental fact 
that HgTiOs crystals are of light yellow color.- 

It is well known that the DFT always underestimates 
the energy band gap. The GW approximationiZiii based 
on the many-body perturbation theory is an approach 
that enables to obtain Eg values in good agreement with 
the experiment. The calculations carried out in this 
work in this approximation gave the band gap energy 
Ef^ w 2.43 eV for the R3c phase of HgTiOg. This 
value agrees with the color of the samples much better 
than the Eg values calculated in the LDA (1.49 eV, this 
work) and GGA (1.6 eV, Ref. [lOl) approximations. The 
optical dielectric constant of HgTiOa which takes into 
account the local field effects is eoo — 9.29. 

To interpret experimental data of infrared (IR) and 
Raman studies, the phonon frequencies calculated at the 
r point for different phases of HgTiOa may be useful. 
The data for R3c and R3 phases at ambient pressure as 
well as for the Pbnm phase at 147 kbar are presented in 
Table FVl In the low-temperature ferroelectric R3c phase, 
Ai and E modes are active both in IR and Raman spec- 
tra. In the i?3 phase, A^ and E^ modes are IR active. 



TABLE V. Calculated frequencies of optical phonons at the F point of the Brillouin zone for HgTiOs phases with R3c, R3, 
and Pbnm space groups (the latter — at P — 147 kbar). 



R3c structure 



Pbnm structure 



Mode 



Mode 



Mode 



Mode 



Mode 



Ai 



Aa 



78 
181 
379 
476 
62 
347 
355 
417 
753 



E 



R3 structure 



73 
204 
306 
440 
651 
133 
348 
477 
647 



Ea 



Eu 



81 
121 
139 
165 
274 
312 
443 
495 
515 



94 
189 
311 
442 
565 
148 
256 
371 
452 



Ao 



All 



B 



Is 



113 
144 
277 
417 
462 
559 
65 
74 
108 
141 
303 
375 
498 
539 
80 
103 
139 
350 



B 



l9 



B; 



2g 



B: 



3g 



Bill 



439 


Biu 


530 


502 


B2u 


38 


782 




86 


104 




141 


266 




190 


450 




345 


542 




356 


819 




431 


118 




496 


226 




524 


353 


B3u 


58 


539 




114 


732 




165 


38 




246 


84 




301 


134 




380 


243 




403 


387 




448 


475 




547 



whereas Ag and Eg modes are Raman active. In the high- 
pressure Pbnm phase (Sec. lIIID|) . the i?i„, i?2u, and i?3„ 
modes are IR active, whereas Ag, Big, B2g, and B^g are 
Raman active. 

The elastic moduh tensor in the R3c phase is pre- 
sented by seven independent components: Cn = C22 = 
348.0 GPa, C33 = 260.3 GPa, C12 = 178.5 GPa, C13 = 
C23 = 149.9 GPa, C44 = C55 = 76.3 GPa, Cge = 
84.8 GPa, and C14 = -C24 = Cgg = 18.3 GPa. The bulk 
modulus calculated from these data is B = 205.8 GPa, 
which is slightly higher than the 178 GPa value obtained 
in Ref. lIO without taking into account the relaxation of 
internal degrees of freedom. 

In the i?3c phase, seven independent components of 
the elastic tensor are: Cn — C22 — 293.9 GPa, C33 = 
224.4 GPa, C12 = 159.7 GPa, C13 = C23 = 117.1 GPa, 
C44 = C55 = 53.7 GPa, Cee = 67.1 GPa, and C14 = 
-C24 = C56 = 0.82 GPa. The bulk modulus in this 
phase is B = 171.3 GPa. 



D. Pressure-induced phase transition 

To discuss the experimental data on the influence of 
hydrostatic pressure on the structure of HgTi03)^i^ cal- 
culations of the pressure effect on the properties of these 
crystals were performed. At pressure P ^ 0, the ther- 
modynamically stable phase is the phase which has the 
lowest enthalpy H = Etot + PV, not the lowest total en- 



ergy -Etot- To compare the phases with different number 
of atoms in the unit cell, we use the specific energy and 
the specific volume defined per one formula unit. The 
calculations show that with increasing pressure, the con- 
tribution of the PV term gives about 95% of the change 
in H in our crystals, and so at high pressures the phase 
with a minimum specific volume becomes more stable. 
As follows from Table [III at P = the Pbnm phase has 
the minimum specific volume. Other phases, arranged in 
order of increasing their specific volumes, form the fol- 
lowing sequence: Imma, lA/mcm, PA/mbm, RSc, RSc. 
The ilmenite structure i?3, which has the minimum total 
energy at P = 0, is characterized by the largest specific 
volume. This enables to expect that as the pressure is 
increased, the ground-state structure will change in the 
following sequence: R3 — ;• R3c -^ Pbnm. Moreover, 
with increasing pressure, the suppression of ferroelectric- 
ity should be observed (the R3c -^ R3c phase transition). 

The enthalpy differences between the phases under 
consideration and the R3c phase as a function of pres- 
sure are plotted in Fig. [21 It is seen that HgTiOs should 
undergo the R3 — >■ i?3c phase transition at P = 38 kbar 
and the R3c -^ Pbnm one at 141 kbar. As both phase 
transitions are accompanied by abrupt changes in the 
volume of the unit cell at the transition pressure (6% 
and 0.71%, respectively), they should be of the first or- 
der. The Imma phase, whose enthalpy at P = is 
lower than that of the Pbnm phase, at higher pressure 
becomes thermodynamically less stable and so can be ex- 
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FIG. 2. The enthalpy differences between _R3, -R3c, Imma, 
and Pbnm phases and the R2>c phase of HgTiOa as a function 
of hydrostatic pressure. 
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FIG. 3. The square of the A2u mode frequency in the _R3c 
phase of HgTiOs as a function of hydrostatic pressure. 



eluded from eonsideration. Similar behavior, when sta- 
ble i?3 phase transformed under high pressure and tem- 
perature to the pressure-stabilized Pbnm phase, from 
which the metastable R3c phase appeared after releas- 
ing the pressure, was observed in MnTiOs?^ FeTiOs^^ 
and ZnGeOajS Ferroelectric LiTaOs also exhibited the 
R3c — > Pbnm phase transition at high pressures,^- 

To determine the pressure of the i?3c — > R3c phase 
transition accurately, the square of the A2u phonon fre- 
quency in the i?3c phase was plotted as a function of 
hydrostatic pressure (Fig.[21). It turned out that this de- 
pendence is non-monotonic and the pressure of the phase 
transition is ^152 kbar. This means that in the whole 
pressure range before the phase transformation to the 
Pbnm phase at 141 kbar the ground state at T = is 



ferroelectric (the R3c phase). In the Pbnm phase, the 
ferroelectric instability is absent. 

An interesting feature which is seen in Fig. |3] is the 
reentrance of the ferroelectric instability in the _R3c phase 
at pressures above 290 kbar. Similar behavior of the fer- 
roelectric mode frequency with increasing pressure was 
predicted for a number of oxides with the cubic per- 
ovskite structurc^'^ and explained by an increased mixing 
of the Ti 3d and the O 2s orbitals at high pressures. The 
effect observed in HgTiOs needs further investigation, 
but it should be noted that it is observed in a crystal 
in which the oxygen environment around the Ti atom is 
different from that in cubic perovskites. 

The obtained results enable to propose new interpre- 
tation of the phase transition observed in high-pressure 
X-ray experiments.— i^ The specific energy of the cubic 
(PmSm) phase, which was considered in Refs. |a and [3 as 
a high-pressure phase, is about 1 eV higher than the spe- 
cific energy of the R3c phase, and its specific volume is 
higher than that of the R3c phase f Table HI]) . This means 
that the large difference in enthalpies of these phases will 
only increase with increasing pressure. Therefore, the 
PmZm phase should not be considered as a high-pressure 
phase. According to our calculations, the pressure coef- 
ficient of the rhombohedral angle in the i?3c phase is 
da/dP = -h0.00547kbar, and so at P = 20 kbar the 
structure remains strongly distorted, with a relative de- 
crease of interplanar distance of about P/3B w 0.32%. 
This value is several times smaller than the relative de- 
crease of interplanar distance observed at the transition 
pressure. 

However, if one assumes that in Refs. |8| and |9| the 
pressure was measured incorrectly (according to our 
estimates, it was underestimated by 5-7 times), and 
one takes the relative change of the (io24 interplanar 
distance^^ as a measure of pressure, the agreement be- 
tween our calculations and experiment becomes satisfac- 
tory. Indeed, at the i?3c — ^ Pbnm phase-transition pres- 
sure (141 kbar) the calculated decrease of the do24 inter- 
planar distance (compared to P = 0) is 2.0%, whereas 
in experiment it is 2.3%. The calculated drop in mean 
interplanar distance^ at the phase transition (0.054%) is 
also close to that observed in the experiment (~0.05%). 

The lattice parameters and atomic coordinates for the 
Pbnm structure at 141 kbar are given in Table IIIIl The 
calculated diffraction patterns for the P3c and Pbnm 
phases at 141 kbar are shown in Fig. |4l The pattern 
for the Pbnm phase is really close to that observed in 
the high-pressure experiment ;^i2i When transforming to 
the high-pressure phase, the (012) line becomes broader 
because in the orthorhombic phase a pair of close lines 
with indices (002) and (110) appears. The (021) line 
characteristic for the orthorhombic phase is clearly seen 
in experimental diffraction patterns obtained during re- 
leasing the pressure!^ The most serious disagreement 
between our calculation and experiment consists in the 
absence of (111) line of the orthorhombic phase in the 
diffraction patterns at high pressure. Possibly, this is 
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first-principles calculations of the thermodynamic stabil- 
ity of mercury titanate were performed. 

To check the thermodynamic stability of HgTiOs , the 
enthalpy of the i?3c phase was compared with that of 
the mixture of starting components, orthorhombic HgO 
and rutile Ti02- The calculations showed that at P = 
the enthalpy of mercury titanate is 150 meV (per for- 
mula unit) higher than the sum of enthalpies of HgO 
and Ti02. This means that at P = mercury titanate 
is thermodynamically unstable against its decomposition 
into starting components. However, because the specific 
volume of the HgTiOa unit cell is significantly lower than 
the sum of specific volumes of HgO and Ti02, the sta- 
bility of HgTiOs increases with increasing pressure. For 
example, at 58.8 kbar the enthalpy of HgTiOs is 75 meV 
lower than the sum of enthalpies of HgO and Ti02. This 
explains why the synthesis of mercury titanate is possible 
only at high pressures. 



FIG. 4. Calculated diffraction patterns for _R3c and Pbnm 
phases of HgTiOs at P = 141 kbar (for Cu A'a-radiation). 



IV. CONCLUSIONS 



due to the incompleteness of structural transformation. 
We think that new high-pressure experiments are needed 
to check the proposed interpretation of the high-pressure 
phase transition in HgTiOa. 



E. Thermodynamic stability of HgTiOs 



As was mentioned in Sec. IHICi the inconsistency be- 
tween the observation of SGH signal and the absence of 
dielectric hysteresis loops in HgTiOa at 300 K can be 
explained by easiness of the defect formation. Indeed, 
according to Ref. 0, the samples darkened when exposed 
to light. The absence of a sharp peak on the temperature 
dependence of dielectric constant at the Curie tempera- 
ture can also be explained by the existence of defects. 
To clarify why the defect formation in HgTiOs is so easy. 



First-principles calculations within the density func- 
tional theory have revealed that the P3c structure ex- 
perimentally observed in HgTiOa is metastable at am- 
bient pressure. With increasing hydrostatic pressure, 
the ground-state structure changes following the R3 ^■ 
P3c -T- Pbnm sequence, and so a new interpretation of 
the phase transition observed at high pressure is pro- 
posed. It is shown that the appearance of ferroelectric- 
ity in HgTiOa at P = is associated with an unstable 
phonon mode. Optical and elastic properties of differ- 
ent phases of mercury titanate are calculated. The band 
gap obtained in the GW approximation {Eg = 2.43 cV) 
agrees with the experimental data better than the value 
obtained in the LDA approximation (1.49 eV). Analysis 
of the thermodynamic stability explains why the synthe- 
sis of mercury titanate is possible only at high pressures. 

The calculations presented in this work have been per- 
formed on the laboratory computer cluster (16 cores). 
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Absolute values of do24 at P = shown in Fig. 3 of Ref. |8| 
disagree with the lattice parameters given in this paper for 
the same pressure (the deviation is about 5%). 
At the i?3c — > Pbnm phase transition the (012) peak splits 
into two components with (110) and (002) indices. 



In the -R3 phase, the dielectric constant is considerably 



